Dilepton production in heavy-ion collisions 
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' The in-medium spectral functions of p and lo mesons and the broadening of nucleon res- 

^ i' onances at finite baryon density are calculated self-consistently by combining a resonance 

, dominance model for the vector meson production with an extended vector meson domi- 

nance model. The influence of the in-medium modifications of the vector meson properties 
on the dilepton spectrum in heavy-ion collisions is investigated. The dilepton spectrum is 
generated for the C-l-C reaction at 2. OA GeV and compared with recent HADES Collabo- 
ration data. The collision dynamics is then described by the Tiibingen relativistic quantum 
molecular dynamics transport model. We find that an iterative calculation of the vector 
y ' meson spectral functions that takes into account the broadening of the nucleon resonances 

' due to their increased in-medium decay branchings is convergent and provides a reasonable 

description of the experimental data in the mass region 0.45 < Af < 0.75 GeV. On the other 
side, the theoretical calculations slightly underestimate the region < M < 0.4 GeV. 
^ ' Popular in-medium scenarios such as a schematic collisional broadening and dropping vector 

, mesons masses are discussed as well. 
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PACS numbers: 12.40Vv,25.75Dw 



I. INTRODUCTION 

It is a well-established fact that hadrons change their properties in a dense and excited nuclear 
medium. Such changes are reflected in mass shifts and/or in the development of complex spectral 
^ ' properties. A typical example is the nucleon that suffers a substantial mass shift at finite density, 

^ ■ see, e.g., Ref. but maintains its good quasiparticle properties. Besides a shift of the pole mass, 

resonances, both nucleonic and mesonic, have the tendency to be broadened and to develop spectral 
distributions that may even lead to a loss of good quasiparticle properties. For example, total 



photoabsorption cross sections on heavy nuclei [2, l3( provide evidence for a substantial collisional 
broadening or melting of nucleon resonances inside the medium Q]. 

To study the medium modifications of hadrons is of particular interest, since it not only pro- 
vides insight into the properties of the strongly interacting hadronic many-body systems but also 
allows conclusions to be drawn on QCD "observables" that characterize the medium. A prominent 
example is the scalar quark condensate {qq) which determines the chiral symmetry breaking scale 
of QCD in the nonperturbative sector. 

Heavy-ion reactions present therefore a unique opportunity for the study of nuclear or hadronic 
matter under extreme conditions, i.e., at supranormal densities and high temperatures. Photo- 
or hadron-induced reactions on the nucleus provide complementary information on cold matter at 
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moderate densities. 

The light vector mesons p and lo are both of particular interest, because their decay into 
dileptons allows one to probe the electromagnetic response of the medium. For this purpose, 
electromagnetic probes such as dilepton pairs have proven to be most efficient, since they leave the 
medium essentially undistorted by final-state interactions. In heavy- ion reactions, they provide a 
clear view of the effective degrees of freedom at high baryon density and temperature. 

Theoretically, an abundance of models can predict the changes of vector meson masses and 
widths in high density, high temperature nuclear matter, and they can be roughly divided in four 
different classes: Brown- Rho scaling Q, models based on QCD sum rule s 
relations - - - -- -~ 
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13], and effective hadronic models [U 



], Q, [13, H, [li [3, [lis 
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The first approaches to the description of the in-medium vector mesons were based on effective field 



theories (EFTs) 2j] and the Nambu-Jona-Lasinio (NJL) model [25|. In some aspects, the various 
approaches come to qualitatively similar conclusions; however, the overall situation is still unclear. 
While Brown-Rho scaling, at least in its naive form, predicts a common downward mass shift 
of the vector mesons where the quasiparticle properties are essentially maintained, the hadronic 
models come to different conclusions. Concerning the p meson, these sets of models predict in 
general a significant broadening of the p and the development of complex structures in the spectral 
functions, e.g., the appearance of additional peaks caused by the coupling to nucleon resonances. 
In some cases, this occurs in line with a slight shift of the quasiparticle peak which corresponds to 
an additional mass shift [lO|. Concerning the lo meson, the situation is even less clear. Early QCD 



sum rules calculations predicted even a repulsive mass shift while in Refs. [3, l9| the strong 
dependence of the oj properties on the higher order unknown quark condensates has been pointed 
out, which leaves room for mass shifts in both directions. The hadronic approaches predict in 
common an essential broadening of the w, although they range from a strong downward mass shift 
[lo| to a slight upward mass shift 2C, 23| to an essential repulsive mass shift jl9| . 

However, recent progress from the experimental side allows one, at least partially, to constrain 
the various theoretical models. While the CERES 26|, |27| and HELIOS 2a] dilepton experiments 
at the CERN Super Proton Synchroton (SPS) revealed clear evidence for in-medium effects in 
heavy-ion reactions (Pb-|-Au) through the observed enhancement of the dilepton spectra below the 
the p and to peaks relative to standard hadronic cocktail sources, such a behavior could be explained 
either within a scenario of a dropping p vector meson mass 29| or by the inclusion of in-medium 
spectral functions for the vector mesons [l3]. Thanks to unprecedented resolution, the recent 
NA60 dimuon experiment 30] was able to "measure" the in-medium p spectral function under 
the conditions of ultrarelativistic heavy- ion collisions. NA60 seems to rule out a naive dropping 
mass scenario but supports the picture of modified p-uj spectral functions predicted by hadronic 
many-body theory 311 ]. 

A second set of heavy- ion experiments have been performed at laboratory energies of 1.0^ GeV 
(Ca-|-Ca and C-l-C) by the DLS Collaboration at the LBNL Bevelac 3j, |33]. Also in this case, 
the low mass region of the dilepton spectra is underestimated by present transport calculations, 
in contrast to similar measurements (1.04-4.88 GeV/nucleon) for the p+p and p+d systems. As 
opposed to the ultrarelativistic case, the situation does not improve when the in-medium spectral 
functions or the dropping mass scenarios are taken into account 3J, |35| (the DLS puzzle) . Other 
scenarios such as possible contributions from the quark-gluon plasma or in-medium modifications 
of the r] mass have been excluded as a possible resolution of this puzzle. Decoherence effects 36 ] 
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have proven to be partially successful in explaining the difference between the DLS data and 
the theoretical predictions. However, in this energy regime, which probes the high density, low 
temperature phase, the situation is going to be improved significantly with the already existing 
and forthcoming measurements of the HADES Collaboration at GSI 371. l38l. l39l|. Coniplemented 
are the heavy-ion experiments by 7-nucleus reactions. The CB-TAPS experiment [40[], which 
focused exclusively on the uj meson and reported an enhanced strength below the oj peak, reports a 
broadening of the uj observed in 7-nucleus reactions. Also, the dilepton mass spectrum measured at 
Japan's National Laboratory for High Energy Physics (KEK) in p + A reactions at a beam energy 
of 12 GeV 
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42( 1 revealed an excess of the dileptons below the p-meson peak over known sources. 
However, these data could not be explained within the standard dropping mass scenario and/or 



assuming a significant collision broadening of the vector mesons 43 ] . An enhanced bremsstrahlung 



44, 



45l |. will most likely not help 



contribution, which is presently under debate at low energies 
explain the high energy KEK data. 

A major difficulty in the interpretation of heavy-ion collision experimental data lies in the fact 
that the gap between observables and theoretically predicted in-medium properties of hadrons has 
to be filled by transport models. Transport models account for complicated reaction dynamics and 
provide the link between theory and experiment. A drawback of such a procedure is dependence on 
the phenomenology and an extended set of input parameters entering the models. Usually, a signif- 
icant fraction of such parameters can be neither constrained by data nor based on well-established 
theoretical approaches. This certainly diminishes the possibility of testing experimentally theory 
and drawing physical conclusions from the experimental data. 

In the present work, we remove some model uncertainties by applying a unified description of 
vector meson production, vacuum decays, and in-medium properties of vector mesons. For this 
purpose, we use a resonance dominance model for nucleon-nucleon scattering in combination with 
an extended vector meson dominance (eVMD) model. Nucleon resonance dominance (NRD) is 
an effective principle which assumes that vector meson production runs over the excitation of 
nucleon resonances 20, |46|, l47|, |48|]. On the other hand, eVMD introduces radially excited p and 



u! mesons 



49| in the RN'y transition form factors 46(] in order to fulfill the quark counting rules 
as a strict consequence of QCD [50]. This allows the kinematically complete, gauge invariant, 
fully relativistic, and unified description of the nucleon resonance transition amplitudes R NV 
{V = uj,p), R A^7, 7*A^ — > R (electro-production), and R Ne~^e~ with arbitrary spin and 
parity in terms of the magnetic, electric, and Coulomb transition form factors. The eVMD model 
solves a long-standing problem of VMD which underestimates the p-meson branchings of nucleon 
resonances when the normalization to the photon branchings is performed. The parameters of 
eVMD are fixed by fitting to photoproducion and electroproduction experimental data, by using 
results of the vrA^ multichannel partial-wave analysis, and, when the experimental data are not 
available, by using predictions of the quark models [46|]. 

Once the model parameters are fixed, one obtains a unified (and parameter free) description 
of quite a broad range of physical processes including vector meson decays, nucleon resonance 
decays to vector mesons and dileptons, and vector meson and dilepton production in elementary 
and heavy-ion reactions. The NRD-l-eVMD model has successfully been applied earlier to vector 



meson (tu and (p) production in elementary {p + p) reactions 5l|, [54] and dilepton production in 
elementary p + p and p + d reactions j47^. In Ref. 53 1 it has been further demonstrated that this 
model is qualitatively able to explain the uj and (p angular distributions in p+p reactions |54l. l55l|. 
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Embedded within the framework of the Tubingen relativistic quantum molecular dynamics 
(RQMD) transport model 56|, |57|, the NRD+eVMD model has been applied to heavy-ion reac- 



tions without introducing new parameters 36|, [S^]- The comparison with the dilepton data from 



DLS 32, ISSJ and HADES [38t | collaborations revealed clear evidence for the in-medium effects 
required, in particular, to suppress excessive dilepton production from the w-meson decays. In 
Refs. 36|, l58|] the collisional broadening and dropping the vector meson masses have been ana- 



lyzed phenomenologically. In the present work, we go beyond the phenomenological analysis by 
calculating the in-medium spectral functions of the p and iv mesons and nucleon resonances using 
the NRD-l-eVMD approach. This allows the first self-consistent theoretical description to be made 
of dilepton spectra based on a unified model for nucleon resonances, vector mesons, and dilepton 
production, and their in-medium modifications. 



II. IN-MEDIUM SPECTRAL FUNCTIONS 
A. Resonance model 

The in- medium properties of hadrons are generally expressed in terms of the self-energy Sy. 
The self-energy determines the spectral function of the quasiparticle in the medium. As long as the 
self-energy shows only a moderate energy dependence, the real part of can be interpreted in 
terms of a mass shift, while the imaginary part generates the in- medium width. To leading order 
in density, the self-energy is determined by the forward scattering length of the hadron with the 
surrounding particles. Since the nucleon and w-nucleon scattering lengths are unknown from the 
experimental side, these quantities have to be determined theoretically. 

In the present work, we apply the resonance model to calculate the forward scattering of vector 
mesons on nucleons. The resonance model is not a field theory in the strict sense where corre- 
sponding Feynman diagrams are evaluated but rather an effective model that has some similarity 
to a field theory based on Feynman diagrams with the intermediate resonances in the s channel of 
vector meson and nucleon scattering. Such an approach was applied in many previous investiga- 



tions of vector mesons properties in the nuclear medium 13j, |20|, |2l|, |59|, |60] . The present approach 
differs with respect to previous investigations by the fact that in the NRD+eVMD model the cor- 
responding couplings of resonances to the nucleon and vector meson are of relativistic form and 
kinematically complete. 

The self-energy of a vector meson V in an isotopically symmetric nuclear medium is deter- 
mined by the invariant VN forward scattering amplitude Ayi\} 

Here V refers either to a p*^ or a cj meson. Due to isosymmetry of the medium, the self-energy 
for mesons is the same as for meson. The forward scattering amplitude Avn is the same 
for proton (A^ = p) and neutron (A^ = n) scattering. The integral in Eq. ([T|) runs over the nucleon 
momenta within the Fermi sphere with Fermi momentum determined by nuclear matter density 

PB 

PB = ^Pf ■ (2) 
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The ampHtude Avn is of Breit-Wigner form for resonance scattering 

. _ (2J/; + 1) Stts ^RNvjs) . s 

^ 2x3 k s-Ml + i^T'-'{s) ' ^ ^ 

In Eq. (I3|) the scattered vector meson has running mass squared and momentum s = {pN+pf' 
is the running mass squared of the baryon resonance i?, and k is the cm. momentum. The width 
rRivy(s) refers to the decay of the baryon resonance R to nucleon and vector meson V with 
fixed mass squared M^. 

The width r^*'(s) refers to the decays of resonance R not modified by the medium, in particular, 
with the vacuum spectral functions for the decay products. This represents the first approximation 
in the calculation of the medium contribution to the total self-energy Sy* = + Sy"* of the 
vector meson V . The vacuum self-energy Sy"* is determined by the corresponding vacuum width 

$js|^^ = -myr*^*(M), = . (4) 

Here T'^p^{M)^ T'^°^{M) are essentially given by the decay widths of the p meson into two pions and 
of the u) meson into three pions, respectively. The two-pion decay width of the p meson is given by 



where k.,^{M,mT^,'m.,^) is the momentum of the pions in the rest frame of the decaying p meson 
having mass M; mp is the physical p meson mass and r*°*(mp) = 150 MeV the on-shell decay 
width. The three-pion decay width of the u) meson can be calculated according to the two-step 



process a; ^ pvr ^ 37r as proposed by Gell-Mann, Sharp, and Wagner [611]. The corresponding 
result can be parametrized in the simple form 

r-(M) = C.(™.) (^5MVe(M^ - (6) 



with m^j the physical w-meson mass, and T'^°^{rn^) = 8.4 MeV the on-shell decay width. 

In the next order, the medium modification of the resonance spectral function including the 
modification of the resonance width due to the modifications of products of the resonance decay 
should be taken into account. 

The width Tjii^y{s) can be expressed by the helicity amplitudes A-i =< 1 — ^\S\^ >, Ai =< 



li|5|i >, 5i =< - > of the R^NV decay [461 



2 



J, 2(yl| 

The calculation of these amplitudes uses the coupling constants of vector mesons to the RN 



transition current. They were obtained in Ref. 46|] by fitting photoproduction and electropro- 



duction amplitudes of baryonic resonances in the eVMD model. The transverse and longitudinal 
self-energies Sy and Sy can be obtained by the following substitutions in Eq. ([7]): 

2,, 9 ,9 ^9^ ,,9 ,9,1 + cos^^ „9 sin^ 6* 

-(^23+^1+51)^(^1+^1) + 25?^—, (8) 

0222 22 Z 2Z 

I (Al +A\+S^)^ 2Sl cos2 e + (Al + A\ ) sin^ 6 , (9) 

3 2 2 2 2 2 2 
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where 9 is the polar angle of vector meson momentum in the cm. system. The polarization 
averaged self-energy Sv reads then 



2Tj'y + 



(10) 



The vector meson spectral function Ay is defined by the off-shell self-energy Sy*(M, |p|) as follows 



Av{M,\p\ 



vr (M2 



m 



V 



)2 + 



(11) 



The helicity amplitudes entering into Eq. 
relativistic approach [46 1 

successfully been applied 

This 



43, 

tering: 



51 



([7]) have been calculated within the same 
and with the same set of baryonic resonances R that has 
to dilepton and vector meson production in p+p collisions 
includes the following set of resonances for pN and ujN scat- 



A^*(1650)i , iV*(1520)| , Af*(1440)i^, iV*(1720)| + , Af*(1680)|^, 



,7 + 



7V*(1535)i" 

A(1620)i~, A(1700)f~, A(1232)|^, A(1905)|"^, A(1950)-2 . 

A straightforward extension of the approach to finite temperature and baryon chemical potential 
would be to integrate the present amplitudes, Eq. ([7j), over hot Fermi distributions. This can easily 
be done and will be a first step toward an application, e.g., at SPS conditions. However, for a 
meaningful determination of spectral functions at SPS conditions, one would have to take into 
account the coupling to not only baryonic but also mesonic excitations (vr, K, . . .) |31]. 

So far, analyticity has not been used in the data analyses to determine the multichannel ttN 
scattering amplitudes 62, l63|, |6j] . The current phenomenological schemes provide resonance masses 
and widths, based on multichannel unitarity and other, less fundamental constraints. The back- 
ground phases entering the dispersion relations are not provided. 

We did not attempt to embed analyticity and restricted our approach to energies s < 4 GeV^, 
where the sum over Breit-Wigner poles gives typically a good approximation for the amplitudes 
65l |. The background is described by t-channel c-meson exchange and the li-channel part of the 
Compton pN scattering diagram evaluated in the Born approximation. 



B. Nonresonant contributions 

Up to now we have not discussed possible nonresonant contributions to the forward vector 
meson-nucleon scattering. The reason is twofold. First, we cannot fix the nonresonant amplitudes 
with the same accuracy as the resonant ones. Second, if we fix them with the available accuracy we 
would find that nonresonant amplitudes approximately cancel in the sum. For example, in the case 
of the p meson, there exist the Compton scattering amplitude which gives a positive contribution 
to the real part of the /?-meson self-energy and the amplitude due to cr-meson exchange which 
gives a negative contribution to it (the latter is of the same origin as the attractive part of the 
A^A^ interaction 66|]). The unknown ppa coupling constant can be extracted from the width of the 
p^ — > 7r"'"7r~7r"'"7r~ decay if one assumes that this decay goes over an intermediate p^a state. 

The two contributions from Compton scattering (5]'-'°™p*) and a exchange (J]"""'^^'^^) are shown 
in Fig. [TJ For the estimate shown in Fig. [T] the corresponding NNp tensor coupling and NNa 
coupling strength were taken from the Bonn one-boson-exchange model f&S] for nucleon-nucleon 
scattering [fNNp = 19.8 (tensor coupling) and gNNa = W ]. The error band for Y,'^-^^'^^ is due to 



7 



the relatively large uncertainty in the four-vr decay of the p meson 

Br(p° 7r+7r"7r+7r") = (1.8 ± 0.9) x 10"^ . 

Nevertheless, from Fig.[Tl one sees that the contributions from Compton scattering and a-exchange 
are of different sign and comparable magnitude. For the mean value of the Br(/9^ — > 7r"'"7r~7r"''7r~) 
branching, they almost cancel completely, and changes of the p-meson spectral function shown 
below are insignificant. 
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FIG. 1: (Color online) Nonresonant contributions to p-meson self-energy from Compton scattering amplitude 
(solid line) and from the amplitude due to exchange by a meson (shaded region). The shaded region 
corresponds to the error in the branching ratio Br{p^ 7r+7r^7r+7r^) = (1.8 ± 0.9) x 10^^. 



To account for nonresonant contributions to the uj spectral function within the present scheme, 
we assume an ujuja coupling three times larger than that for ppa which is motivated by the com- 
parison with the two-pion coupling. The NNlo vector coupling {gNNui = 15.9) is again taken from 
the Bonn potential 66|. As can be seen in Fig. HI the influence of the nonresonant contributions 
is now more pronounced than in the case of the p meson; however, the qualitative features of the 
spectral distributions are not changed. 

The Breit-Wigner amplitudes decrease as 1/s with increasing s. Such a parametrization ensures 
the change of the resonance phases by tt from low to high energies. The c-meson exchange generates 
the scalar mean field, which is known to be important in the modification of the nucleon masses 
[l|. It plays an important role in our scheme too. The component of the amplitude connected to 
the CT- meson exchange remains constant for s ^ oo. 



C. p-meson spectral function 

In the following, we discuss first the p meson. Figure [2] shows the p spectral function in nuclear 
matter at nuclear saturation density po = 0.16 fm~^. Longitudinal (A^) and transverse {A'^) 
spectral functions are found to be rather similar. This means that unpolarized spectral functions 
can be used in the calculations of dilepton spectra. 
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FIG. 2: (Color online) Longitudinal (L) and transverse (T) p spectral functions in nuclear matter at satu- 
ration density for various momenta p (in GeV). Dashed lines stand for the resonance approximation, solid 
lines represent calculations that also included the nonresonant contributions. The shaded area shows the 
vacuum spectral function. 



We observe a slight upward mass shift of the p and a substantial broadening. At low momenta, 
the spectral functions show a clear two-peak structure which vanishes with increasing vector meson 
momentum. The results shown in Fig. [2] are in qualitative and even quantitative agreement with 
previous calculations based on the resonance model assumption |21)]. Although the various ap- 
proaches are based on different ways to describe the corresponding transition form factors, eVMD 
in the present case, and parameters are partially fixed in different way, this fact demonstrates the 
stability of the essential features predicted by these types of models. 

The emerging two-peak structure can be understood as follows. The value and sign of the self- 
energy depend on the pole positions of the particular resonances. If the vector meson mass 
squared is small, the invariant mass of vector meson plus nucleon is below the pole masses of the 
relevant nucleon resonances. Therefore the real part of the vector meson self-energy is negative. 
This is a typical example for level repulsion (vector meson plus nucleon and nucleon resonance). 
Consequently, the factor {m? — niy — SRSy)^ in the denominator of the vector meson spectral 
function, Eq. (jlip . is small or even equal to zero. Thus the first peak in the spectral function 
emerges at a vector meson mass around 0.5 GeV. The major contribution, which generates the 
first peak, comes from the N*{lb2Q), which is in agreement with the findings reported in Ref. 211]. 

If the vector meson mass squared lies in the vicinity of its vacuum value niy, the invariant mass 
of vector meson plus nucleon lies above the pole masses of the relevant nucleon resonances and 
the real part of the vector meson self-energy is positive. Thus we obtain the second peak in the 
spectral function at a vector meson mass slightly above my. 

At high vector meson momenta, the invariant mass of the vector meson plus nucleon is always 
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above the pole masses of the relevant nucleon resonances. As a result, the spectral function has 
only one single peak slightly above my- 
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FIG. 3: (Color online) Unpolarized p meson spectral function at rest in nuclear matter at saturation density 
and twice saturation density. The shaded area displays the vacuum spectral function. 



Figure [3] displays finally the dependence of the p- meson spectral function on nuclear density. It 
shows the unpolarized /o-meson spectral function at rest at pQ and at 2po nuclear density. With 
increasing density, we observe a further shift of strength away from the original pole mass; i.e., the 
first branch in spectral distribution is slightly enhanced and even shifted to lower masses, while 
the second peak is slightly shifted upward at 2pQ compared to pQ and also additionally broadened. 

In this context, it should be noted that the resonance model predictions stand in contrast to the 
EFT coupled-channel calculations of Ref. 19|], which predict no significant medium dependence 



of the p, concerning neither a mass shift nor a broadening. The reason that in the approach of 
Ref. [l9t] much less strength is shifted to lower masses lies mainly in the much weaker coupling to 
the TV* (1520) found in Ref. 0. For this resonance, the value of Tnp ~ 2 MeV fiol] has to be 
compared with F^rp ~ 25 MeV from Refs. 2l|, |46|. The latter value, however, agrees with that of 
the PDG 



671 1 the Manley and Saleski analysis [62 [. 



D. w-meson spectral function 



For the uj meson, we observe a behavior that is principally similar to that of the p meson (see 
Fig. H]). Transverse and longitudinal spectral functions are similar. In both cases, the u pole mass 
is slightly shifted upward, and the uj is substantially broadened around its quasiparticle pole. At 
po we obtain an in-medium u width of 300 MeV. 
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As in the case of the p, the couphng to low lying resonances leads to the appearance of a first 
peak in the spectral function which lies around 0.5-0.55 GeV. With increasing momentum, this peak 
is washed out and disappears finally. However, in the case of the uj, the influence of nonresonant 
contributions is found to be much stronger than for the p. The nonresonant contributions tend 
to increase the repulsive mass shift of the to pole, and they strongly suppress the first peak in the 
spectral function. 

This first branch in the spectral distribution is mainly generated by the A'^*(1535) resonance. 
As discussed in detail in Refs. [4a, [511], within the NRD+eVMD model a strong A^*(1535)A^a; 
coupling is implied by the available electroproduction and photoproduction data. However, the 
Nuj decay of this resonance has not been measured directly, and therefore input from quark model 
predictions had to be used to fix the entire set of eVMD model parameters. Nevertheless, within 
such a procedure, a strong A^*(1535)A^ti; coupling seems practically unavoidable. In pp — > ppco 
production, the large N*{1535)Nuj decay mode leads to substantial contributions in a kinematic 
regime where the to is far off-shell, i.e., at small invariant masses. This is reflected in an enhancement 
in the cross section around threshold j5l|. Existing data 55, 68, 69 1, however, do not rule out such 
a behavior. A closer inspection of the experimentally observed background contributions may 
provide important experimental information concerning this question. 
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FIG. 4: (Color online) Same as Fig. [21 but for uj spectral functions. 



The nuclear matter density dependence of the cj-meson spectral function is shown in Fig. \E[ 
Again, the figure shows the unpolarized spectral function at rest at po and at 2po nuclear density. 
As for the p meson, we observe a shift of the second peak which belongs to the original uj pole 
toward higher masses with an increase in density, while the first peak is slightly shifted to lower 
masses. Moreover, the height of the second peak is suppressed by about a factor of 2. 

Comparing this with other works, we should mentioned that in the pure resonance model 
approach of Ref. 20|, no such additional peak was observed. The w-meson spectral functions 
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obtained within the coupled-channel approach of Ref. and within the coupled-channel K 

matrix of Ref. [23] have qualitative similarity with those from the present approach. All approaches 
come practically to the same conclusions: an upward mass shift, a broadening of the w, and the 
appearance of an additional branch in the u: spectral function. This branch appears at the same 
position and is in both cases generated by the A^*(1535). However, in all approaches, the uj survives 
as a quasiparticle, at least at moderate densities up to po; i-e., there the spectral function is still 
dominated by the main branch corresponding to the original u) pole. The predictions for the density 
dependence of the spectral function are similar on a qualitative level; i.e., when going from one to 
two times nuclear density, the suppression of the branch corresponding to the u: pole is of similar 
size. 

However, on a quantitative level, the models come to different conclusions. While the broadening 
of the (J is similar in Refs.0 and li, the mass shift is much larger in Ref. [l9| (Am^j ~ 46 MeV 
at po) than in Ref. [23|] (Am^^ ~ 10 MeV at po). In the present case, the in-medium modifications 



of the u meson are even more pronounced than in Refs. 19|, l23|]; i.e., the broadening and the 
upward mass shifts are larger (Attt-o; ~ 75 MeV at po)- 

A comparison with predictions from QCD sum rules 0| turns out to be difficult because the 
u) properties depend strongly on higher order condensates. Sum rules leave space for upward and 
downward mass shifts, and the parameters related to the higher order terms in the operator product 
expansion have to finally be fixed from experiments 0]. Moreover, these approaches assume that 
the oj maintains its quasiparticle properties. However, due to the distinct two-peak structure of the 
present spectral distributions, it is not possible to assign a common mass shift to an uo quasiparticle 
pole. 
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E. In-medium resonances: Role of self-consistency 

As the next step, we took into account the changes induced by the in-medium vector mesons 
on the total width of the nucleon resonances. This leads to a self-consistent determination of the 
self-energies of the vector mesons in nuclear matter. 

The results shown in the previous section correspond to the first iteration, if considered in 
the context of a self-consistent calculation. In the second iteration, the in-medium widths of the 
nucleon resonances F^j are determined by insertion of the in-medium spectral functions of the 
vector mesons resulting from the first iteration. Because the latter depend on the momentum of 
vector meson with respect to nuclear medium p, the in-medium widths of the nucleon resonances 
rjj will depend on the resonance momentum \pr\, that is, 

r*j,{s,\pR\)=T'^\s)+J2 I TRMv{s,M)^Av{MM)dM^'^ , (12) 

V 

where A^y refers to the modification of vector meson spectral functions with respect to the 
vacuum ones. |p|, among other things, depends on \pr\ and on the orientation of the decay 
products momenta with respect to the direction of the resonance momentum. 

Doing so, nucleon-resonance scattering terms leading to the broadening of the resonances 



are produced 211]. The vector meson self-energies are then calculated from Eqs. ([T])-(l3]) using 
r|j(s, IphI) instead of r5^'(s). The procedure is repeated until convergence. We find that the 
convergence is obtained after the third iteration. 

As a side result of our self-consistent calculation, we find that the widths of the nucleon reso- 
nances are enhanced in medium because the vector meson spectral functions show a significative 
spectral strength at small invariant masses. A similar outcome emerged from the analysis per- 



formed in Ref. 70(]. 



The resulting unpolarized vector meson spectral functions are shown in Figs. [6] and [7] for the p 
and OJ mesons, respectively. They refer to saturation density. We observe that the self-consistent 
calculation leads predominantly to a reduction of the lower mass peak. This result qualitatively 



agrees with the findings of Ref. [2]| , which investigated the role of a self-consistent iteration scheme 



on the p-meson spectral function. 



F. Experimental situation 



Experimental constraints on the in-medium a; spectral function can presently be derived from 
the CB-TAPS (7-1-^) experiments [40], thep+A measurements at KEK [4_l|,|42], and the heavy-ion 
dilepton experiments. The dilepton measurements of the DLS Collaboration in C-l-C and Ca-|-Ca 
at \A GeV 84] suffer from too low mass resolution in the vicinity of the oj peak in order to make 



precise statements on the oj in-medium width. However, there is no doubt that the explanation of 
the DLS data requires a substantial broadening of the 00 spectral function. The analysis of Ref. 
3a] showed that the DLS data are compatible with a rather large oj width; i.e., rj;^'* ~ 150-300 
MeV. The first data from HADES 33] will be analyzed in the next section. 



As discussed in Ref. 



431 ]. the interpretation of the p+C and p+Cn KEK dilepton data [4ll. |42| 



suffers from the high initial proton kinetic energy of 12 GeV. This means that vector mesons are 
produced with high momenta {p > 1 GeV) and in particular the oj decays at low nuclear densities 
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FIG. 6: (Color online) Unpolarized spectral functions of the p meson in nuclear matter at saturation density 
for various momenta p (in GeV). The broadening of the nucleon resonance widths induced by the in- medium 
spectral properties of the vector mesons is taken into account and a self-consistent calculation is performed. 
The shaded area shows the vacuum spectral function. 
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FIG. 7: (Color online) Same as Fig. [6l but for the uj meson. 



or even outside the target nucleus. Nevertheless, the KEK study observed a substantial difference 
of the dilepton spectrum with respect to the standard sources below the p/uj peak. On a qualitative 
level, these data support a picture as predicted by the resonance model, i.e., a substantial shift of 
spectral strength to smaller masses. 

However, the resonance model and coupled-channel predictions for the uj spectral function 
contradict, at least partially, the results of a recent photoproduction experiment carried out by 
the CBELSA/TAPS Collaboration 40(]. The experiment indicates a downward mass shift of about 
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Arrii^ ~ —80 MeV at po, whereas the coupled-channel [19|, |23| and resonance model calculations 
predict a more or less pronounced repulsive shift ranging from +10 to +80 MeV. The collisional 
broadening extracted by CBELS A/TAPS is moderate; i.e., AF^; 50 MeV. Since this experiment 
has been carried out at beam energies of 0.64-2.53 GeV, a similar argument as for the KEK 
experiment applies, at least for the high energies, namely, that one has to carefully account for the 
energy dependence of the cv spectral function and to follow the paths of the oj decays. 

Recent j+A measurements from the CLAS Collaboration [7l[ at the Thomas Jefferson National 
Accelerator Facility, carried out at photon energies = 0.6-3.8 GeV on light (C) and heavier 
(Fe) targets, find no signatures for a vector meson mass shift for the p and uj mesons but indicate 
a collisional broadening of the p by 50-70 MeV. 

In terms of a simple BR scaling interpretation my = my(l — apB/po) which should hold in 
common for both the p and the co meson, the best fits to the various experiments yield at present 
a divergent picture: \a\ ~ 0.13 (CBELSA/TAPS \^), \a\ ~ 0.092 ± 0.002 (KEK-PS E325 0), 
and |a| ~ 0.02 ± 0.02 (CLAS [Til). 

This apparent contradiction implies already that the mass shift scenario a la Brown-Rho is 
too simple, and, consistent with the NA60 heavy-ion data 30], the vector mesons develop more 
complex spectral properties. 



G. Realization within the transport approach 

A first attempt to introduce in-medium spectral functions of the vector mesons in a transport 



description for intermediate energy heavy-ion collisions was performed in Ref. [3J]. The in-medium 
dilepton rate was thereby expressed in terms of the p meson in-medium spectral function. The 
proportionality was achieved by including the medium effects at the level of production channels. 
However, to avoid double counting, this required switching off the explicit p meson production 
channels included in the self-energy calculations, which have been implicitly accounted for in terms 
of the in-medium spectral function. In particular, the decays of the p mesons produced in baryon- 
baryon collisions and meson-baryon interactions were not included explicitly. 

In the language of a resonance model, this would mean, e.g., that since V + N ^ R is a (dom- 
inant!) contribution to the self-energy, the R NV Ne^e~ decay should not be included 
explicitly, if one would operate as in Ref. [34]. However, nucleon resonances are important dynami- 
cal degrees of freedom of a transport approach. In particular, in heavy-ion collisions at intermediate 
energies, the medium is dominated by nucleons and nucleon resonances, and to neglect an explicit 
(dynamical) treatment of the latter in the determination of vector meson production to obtain a 
direct proportionality to the vector meson spectral function is a questionable procedure. Moreover, 



it was already pointed out in Ref la] that such a treatment [34| might not be realistic. 

Here we adopt an alternative approach to extract information on the in-medium vector meson 
spectral functions from dilepton emission in heavy-ion collisions. The idea is to restrict propagation 
and mutual interactions to the dynamical degrees of freedom within the transport approach. Vector 
mesons and their interactions are treated as perturbative degrees of freedom within the transport 
description. This gives us the possibility of including the corresponding vector meson in-medium 
modifications on a microscopic level. The philosophy behind this approach is similar to that of the 



approach pursued in Ref. la, l31( , which derived local dilepton emission rates from the decay rates 



of the corresponding sources within the framework of an expanding fireball model. 
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To be more precise, within the NRD model the vector meson production channels are nucleon 
resonances which are treated as explicit, i.e., as dynamical degrees of freedom in the transport 
code up their decay. Dilepton emission takes place via resonance Dalitz decays (eVMD) where 
the vector mesons enter as virtual particles the RN'j* vertices. In the medium, the vector mesons 
entering into the vertex form factors are modified by the vector meson self-energy, determined for 
the conditions, i.e., at the density where the resonance decay takes place. Thus the presence of the 
medium changes the branching ratios for the nucleon resonance Dalitz decays. This modification, 
however, is not directly proportional to the vector meson spectral function but rather to the in- 
medium form factors. In this approach, e.g., vector meson absorption processes are taken into 
account microscopically, although not dynamically, in terms of the imaginary part of the vector 
meson self-energy, which determines the corresponding coUisional broadening. 

By this procedure, one avoids the complicated and yet not fully resolved problem of a consistent 
off-shell propagation of the vector mesons within semiclassical transport models 3^,172]. A draw- 
back is certainly that one loses in this picture information on the dynamical propagation of the 
vector mesons. Consequently, their in-medium properties are determined by the conditions at the 
decay points of the nucleon resonances. Since the vector mesons appear only as intermediate states 
in the resonance decay rates, it is, on the other hand, much easier to include spectral functions 
and to keep quantum effects which are lost in semiclassical approaches, even when off-shell effects 
are taken into account [73j, |74 1 . 

The in-medium spectral properties discussed in the next subsection are determined in a local 
density approximation (LDA). In RQMD, the local baryon density p_B is determined by the sum- 
mation over Gaussian wave packets of all nucleons and resonances (it should not be mixed up with 
the interaction density used in RQMD for the determination of the intranuclear forces 75|). Thus 
the determination of the baryon density does not require local equilibrium as in hydrodynamics, 
but dependences on particular nonequilibrium effects such as phase-space anisotropics 76| or mem- 
ory effects [77]] are neglected. As usually done at intermediate energies, an explicit temperature 
dependence of spectral properties is neglected as well. It is, however, possible to extract local 
temperatures from transport simulations, either by fitting hot Fermi distributions to local momen- 
tum space configurations [78] or by performing thermal model fits to local hadron abundances and 
spectra 79|. Although both procedures are connected with an extremely high numerical effort, it 
is thus possible to include density- and temperature-dependent spectral functions into transport 
simulations. Since temperature effects have been found to dominate at SPS energies [sij], such an 
extension will be necessary for an application of the present model to, e.g., NA60 data. 



H. In-medium dilepton emission rates 



Because of the P invariance of the electromagnetic interaction, resonances with arbitrary spin 
have only three independent helicity amplitudes in the 7*A^ — > R transitions. This means that 
there are three independent scalar functions to fix the vertices. The three scalar functions arising 
from the decomposition of the ■j*N ^ R vertex over the Lorentz vectors and the Dirac matrices 
are functions of the mass squared of the virtual photon and are called covariant form factors. 
In the eVMD model, each of these covariant form factors is expressed in a gauge invariant way 
(see Appendix) as a linear superposition of the contributions from the intermediate vector mesons 
of the p and uj family. In contrast to the naive VMD, in which only the p and uj ground states are 
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taken into account, eVMD includes radial excitations p(1450), /?(1700), etc., which interfere with 
the ground-state p mesons in radiative processes. The corresponding transition form factors are 
given by 46 1 



i^f^(M^) = EMf , (13) 

i 

where /c = 1, . . . , 3 stands for each of the form factors, (zt) denotes states of normal and abnormal 
parity, respectively, and the sum is over the intermediate mesons. The A resonance form factors 
have contributions from only the p-meson family, whereas the nucleon resonances receive contri- 
butions from the p and uj mesons. For a resonance of spin J = I + 1/2, the total number of vector 
mesons is / + 3. The amplitude 

Mt} = hif , ."^l (14) 

represents the contribution of the ith vector meson to the form factor of type k. The residues 
^ki contain the free parameters of the model. They are constrained by the requirement that the 



asymptotic expression of the form factors is consistent with the quark counting rules 80(] . For each 
form factor, the quark counting rules reduce the number of free parameters from / + 3 to 2 for k = 1 
and to 1 for A; = 2, 3. The remaining parameters are fixed by fitting the available photoproduction 
and electroproduction data and using results of the multichannel partial-wave analysis of the ttN 
scattering. Where experimental data are not available, predictions of the nonrelativistic quark 
models are used as an input. 

The T{R ^ ^7*) decay width can be written in terms of three transition form factors (magnetic, 
electric, and Coulomb) for a resonance with spin J > 1/2 and two for J = 1/2. The matrix elements 
connecting the former with the covariant form factors are explicitly listed in Ref. |46| . 

In this representation, the insertion of the in-medium properties of the uj and p vector mesons 
is straightforward. In the medium, the transition amplitudes -M^j^^ {i = p, to, . . .) are directly 
modified by the in-medium self-energy and read 

We include the self-energy contributions for the ground-state p and lo mesons in the transition. 
For the excited states . • - j the self-energies are unknown, and thus we keep for these states 

their vacuum properties. 



As in Ref. 5a], we also consider scenarios in which the self-energy is based on different model 
assumptions, namely, a simple Brown- Rho (BR) or Hatsuda-Lee scaling of the vector meson masses 

0] and a collisional broadening of the vector meson widths. In the latter case, the self-energies 
are given by 



= -mv(ri,°)(M)+rf'(/,B,M) 
5iS*/ = . (16) 



In this context, we want to stress that in Eqs. (jl6p . the energy dependence due to the two- 
or three-pion decay of the vector meson is kept in the vacuum contribution to the total width, 
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while the coUisional broadening due to the interaction with the surrounding nucleons is absorbed 
into a density- and energy-dependent part. The issue of the energy dependence of the coUisional 
width will be discussed in detail in the next section. The BR scaling is introduced through the 



replacement my — > niy = mi/(l — ct^), as done, e.g., in Ref. [8l|]. In particular, in this case, one 
has 

3f?SV°* = (mv- a^) -m\ . (17) 
V Po / 

As usual, the mass shift entering into the real part can be adjusted by the parameter a. As in the 
case of full spectral functions, the self-energy components enter into the amplitudes (jlSp . In this 
context, it is important to note that the modification of the amplitudes (jlSp leads to a coherent 
summation of the p and uj spectral functions in the transition form factors ()13p . Doing so, this 
approach goes beyond the standard-even off-shell-transport approach where spectral properties 
are treated at the level of cross sections [3J, |7^, |73|. The latter always leads to an incoherent 



summation of the contributions from different hadrons. 

The self-energy appearing in Eq. (jlSp is a function Sy'*(M, |p|, ps) of the vector meson running 
mass, the modulus of its three-momentum in the nuclear matter rest frame, and the local density 
of the surrounding matter. In the rest frame L* of a resonance R with mass /i, decaying into a 
nucleon and a vector meson of mass M, the modulus |p*| of the momentum of the meson is fixed 
by energy conservation. If is the momentum of the resonance R in the cm. frame L of the 
colliding nuclei and v/j = p/j/ ^ + is its velocity, the vector meson momentum in L is given 
by the Lorentz transformation 

|2 I XT'* I ^, ^* \2 I 



pr = {ir\^r\e* -riRvlY +p¥ , (18) 



where 



pI = |p*|cos^ , (19) 
Pt = |p*|sin0 , (20) 

with 6 being the polar angle of the meson in L* if one chooses the z axis of this frame pointing 
in the direction of v/j. Since |p*| is fixed, in terms of the L frame variables, one has Sy* = 
T,y^ (M, cos 9, pb), and the decay amplitude averaged over the angles reads 

r+i dcos(9) 

T{R^Nj*){p,M,pb) = -^^r{R^ Nj*){p,M, cose, pb) . (21) 

Equation ()2ip is implementable in the framework of the Tubingen RQMD transport code. The 
RQMD code 0, 0, Q has been extended to include all nuclear resonances with masses below 
2 GeV, in total 11 N* and 10 A resonances. A full list with the corresponding masses and decay 



widths to various channels can be found in Tables III and IV of Ref. 36f|. For each resonance, 
RQMD provides the values of the three-momentum components (necessary to perform the Lorentz 
boost), the mass (distributed over a Breit-Wigner) , and the local density of the surrounding matter 
at the decay point. 

Since vector mesons play in the eVMD model the role of intermediate virtual particles, their 
off-shellness is fully taken into account in a consistent manner. 

The model can be applied to dilepton production in heavy-ion reactions. In the energy range 
of a few A GeV, one can identify three main classes of processes that lead to dilepton emission: 



18 



nucleon-nucleon bremsstrahlung, decay of light unflavored mesons, and decay of nucleon and A 
resonances. Dilepton production through the bremsstrahlung mechanism has been studied in 
detail in Ref. 84]. For the energy range of interest in this work, bremsstrahlung contributes 
in a significant way only at small invariant masses to the dilepton spectrum. By far, the dominant 
contributions result from diagrams that involve the excitation of an intermediate A resonance. 
Within the present framework, the inclusion of such contributions would, however, lead to a double 
counting and therefore we omit explicit bremsstrahlung contributions. Recently, the quantitative 



importance of bremsstrahlung contributions has again been discussed in Ref. [4J], however, with 
results contradictory to those in Ref. [82]. 

At incident energies of a few A GeV, the cross sections for meson M = r],!]' , p,uj,(f) production 
are small, and these mesons do not play an important role in the dynamics of heavy-ion collisions. 
Their production can thus be treated perturbatively, in contrast to the case of the pion. The decay 
to a dilepton pair takes place through the emission of a virtual photon. The differential branching 
ratios for the decay of a meson to a final state Xe~^e~ can be written 



with ^ the meson mass and M the dilepton mass. Three types of such decays have been considered: 
direct decays M. — > e'^e" , Dalitz decays M. 7e+e~ and M. — > 7r(7/)e"*"e~, and four-body decays 
M — > 7r7re"'"e~. A comprehensive study of the decay of light mesons to a dilepton pair has been 
performed in Ref. [49| . Assuming a NRD model for the production of p and uj mesons, the remaining 
decay channels that are most important quantitatively for heavy-ion collisions at \A and 2A GeV 
are vr'^ — > ^e^e~ and t] — > '^e^e~ . 

In terms of the branching ratios for the Dalitz decays of the baryon resonances, the cross section 
for e^e~ production from the initial state X' together with the final state NX can be written as 

dM^ -'^Jim^+MY ^ dp^ ^ dM^ ■ 

Here, p is the running mass of the baryon resonance R with the cross section da{s, p)^''^'^^, 
dB{p, M, pb)^^^^^^'^^'^ is the differential branching ratio for the Dalitz decay R Ne^e~ 
through the vector meson V . Thus Eq. (j23p describes baryon- and pion-induced dilepton produc- 
tion; i.e., the initial state can be given by two baryons X' = NN, NR, R'R or it runs through 
pion absorption X' = irN. In the resonance model, both processes are treated on the same foot- 
ing by the decay of intermediate resonances. Medium modifications enter the branching ratio 
dB{p, M, pb)^~*^^~*^^^''~ by affecting the Dalitz decay width dT{p, M, pb)^~'^^^^''^''~ /dM'^ . 



Once the T(R N'j*) is calculated within the eVMD, the factorization prescription 4a] can be 
used to find the dilepton decay rate 

dT{R ^ iVe+e~) = T{R ^ Nj*)MT{j* ^ e+e")^^, (24) 

where 



Mr(7* - e+e-) = ^(M^ + 2ml)Jl- ^ (25) 



3' "'V 

is the decay width of a virtual photon 7* into the dilepton pair with the invariant mass M 
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As discussed in Ref. [21[, the excitation of particle-hole pairs in the meson spectral function 
generates resonance-nucleon scattering terms in the resonance self-energy and thus the in-medium 
broadening of the resonance. We have mentioned that nucleon resonances are dynamically treated 
in the RQMD model, and resonance-nucleon scattering is explicitly performed. Thus, the in- 
medium broadening of nucleon resonances is taken into account in the transport approach dynam- 
ically. No in-medium spectral functions of the vector mesons are therefore included in the total 
width r^^-(^) entering in the branching ratio in Eq. ()23l) . 

The real part of the resonance self-energy 3?S^^ is included in a phenomenological way through 
a mean field in which the resonances are propagated. It is therefore assumed that the nucleon 
resonances feel the same potential as nucleons, which is a standard approximation in present 
transport models and should be improved in the future. For example, the RNV vertex gives rise 
to additional Fock contributions which could be included in future work. 

An observable tightly connected to a correct treatment of the resonance dynamics in heavy-ion 
collision transport calculations is provided by the pion multiplicity. For the mass system under 
consideration, pion multiplicities are reasonably well reproduced by the present description. For 
example, inclusive vr^ cross sections in C+C reactions measured by the KaoS Collaboration 83 1 
can be reproduced by the present description within error bars. This gives, at least on a global 
level, manifest credit to our treatment. 

For the r/, we include r] absorption from the dominating channel r] + N ^ A'^*(1535) explicitly. 
Since chiral perturbation theory predicts practically no modifications of the in-medium rj mass |84l |. 
we do not include a possible rj mass shift. 

To give an impression of the density range relevant for dilepton production in the C+C system. 
Fig. [8] shows the density distribution dN/dpB where the nucleon resonance decays into dilepton 
channels take place. Note that Fig. [8] refers to minimal bias conditions. Triggering on central 
reactions and/or increasing the system size will help to better explore the high density range. 
To have a separate look at p and oj production, we distinguish between A^*, A*, and A(1232) 
resonances. It can be seen that the highest mass resonances, i.e., N* and A*, decay at supranormal 
densities, while a large fraction of the A(1232) decays take place at lower densities, between 0.5/?o 
and Ipo- In all cases, however, maximal densities up to SpQ-Apo are reached. Thus, already the 
small C+C system probes the spectral properties of intermediate vector mesons from A^* and A* 
decays at supranormal densities. 



III. DILEPTON PRODUCTION 

In this section, we provide theoretical calculations of the dilepton emission in heavy-ion collisions 
at intermediate energy. In particular, we address the reaction C+C at 2A GeV for which exper- 
imental data have been already released by the HADES Collaboration. The main purpose is to 
compare calculations that include in-medium effects in a more traditional way, i.e., via Brown-Rho 
scaling of the vector masses and empirical collisional broadening of the decay width, with results 
obtained using p and uj mesons described by the in-medium spectral functions of the previous 
section. New HADES data 85|] will be analyzed elsewhere. 

In the transport calculation, the reaction has been treated as minimal bias collisions with 
respective maximal impact parameter 6max = 6.0 fm. For the nuclear mean field, a soft momentum- 



dependent Skyrme force {K = 200 MeV) is used 75| which provides also a good description of the 
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FIG. 8: (Color online) Impact parameter weighted density distribution of the nucleon resonances taken into 
account in the calculations at their decay point. The distribution of the N* , A*, and A(1232) resonances 
are separately shown. 



subthreshold production in the considered energy range [8g]. To perform the comparison with 
the HADES data, dilepton events originated from the different considered sources were generated in 
the phase space. After smearing over the experimental momentum resolution, the acceptance filter 
function provided by the HADES Collaboration was applied. Events with opening angle Og+e- — 9° 
were rejected, in accordance with the treatment of the experimental data. The spectrum was then 
normalized to the corresponding vr'^ multiplicity. 



1. Vacuum 



We start by addressing the results obtained without any additional medium effects concerning 
the dilepton production. In Fig. O the dilepton spectrum obtained within the vacuum formulation 
of the NRD+eVMD model is compared with the HADES data [38]. The experimental data are 
slightly underestimated in the mass region rriT^ < M < 0.4 GeV and overestimated in the region 
of the vector meson peak. Indeed, the comparison with DLS data had already shown that the 
eVMD model in its pure vacuum formulation fails in describing dilepton production in heavy-ion 
collisions [36] . However, the vacuum calculation is a good reference point for isolating, where 
possible, those sources that dominantly contribute to the spectrum in a certain invariant mass 
region. Once the dominant sources have been individuated, it is interesting to look separately at 
their modifications due to in-medium effects. For this purpose, we also show separately in Fig. [9] 
the contributions to the spectrum of the decays of the pseudoscalar rj and vr'^ mesons and all the N* 
as well as the A resonances. In addition, the A(1232) iVe^e" decay channel is explicitly shown. 
In what follows, we will investigate the modification of the Dalitz decays of the baryon resonances 
due to the introduction of the in-medium properties of the p and lo mesons. Since we introduce 
no in-medium modifications of the vr'' 7e"^e~ and tj^ — > 'ye~^e~ channels, the contribution to 
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the dilepton spectrum from the vr*^ and rj Dahtz decay will remain unchanged in the course of our 
analysis. 

The rj multiplicity M^(47r)[10~^] = 330 for the considered reaction C+C at 2 A GeV is in 
agreement with experimental data from TAPS 87|: M^(47r)[10~^] = 294 it 46. 

Before coming to the discussion of medium effects, the vacuum results should be examined more 
closely. As compared to that reported in Ref. 45|], we find a higher yield around the p/uj peak 
region. This enhancement arises in the present model first because of a strong coupling of the to to 



the A''* (1535) resonance as discussed in detail in Refs. [36|, l5a]. Second, additional enhancement 
results from the implementation of the quark counting rules to the nucleon resonance transition 
form factors. The quark counting rules are known to be a well-founded consequence of QCD 
and, furthermore, are required experimentally to match the photon and p-meson branchings of the 



47|. 



nucleon resonances also 

A precise estimate of the u; contribution is particularly important for extracting the u; meson 
in-medium width: the underestimation of the dilepton yield gives rise to the underestimation of 
the width. As we shall see, the in-medium u peak is strongly suppressed because of the u meson 
broadening. 

The low-mass region is critical for understanding the DLS puzzle. In the present vacuum calcu- 
lation, we obtain a low-mass dilepton yield that is about a factor of 2 smaller than that in Ref. 45l ]. 
Bremsstrahlung cannot explain this deviation, since at 2A GeV it is marginal whatever maximal 
cross sections under debate fiJ] would have been used. At small M, the variances in predictions of 
the present transport models arise from two additional sources, namely, the rj contribution and the 



A Dalitz decay: For the tj decay which dominates the low-mass dilepton yield [36l. I45l. I58l. I88l. 1891] . 
we obtain quite standard values. The main difference lies therefore in the A Dalitz decay. The 
present calculation is close to that in Ref. 88|] and about a factor of 5 lower than the calculation 
of Ref. 45(. This point is crucial, since the whole interpretation of the low-mass dilepton spectra 
depends on this fact. 

The problems on the A(1232) Dalitz decay have occurred already at a kinematic level where 
theoretical calculations of several groups surprisingly disagree with each other pairwise (for a de- 
tailed discussion, see Ref. [90|]). The dilepton decays can be determined from the radiative decays 
by factorization, which means that the A^7* i-^ A amplitudes have to be determined first. In this 
context, it should be noted that from all nucleon resonances, the A^7* ^ A(1232) transition ampli- 
tudes are the best constrained from the experimental point of view (see, e.g.. Fig. 20 in Ref. j46l|). 
The A(1232) Dalitz decay is dominated at M ~ by the magnetic form factor. The normalization 
at M = 0, assuming the dominance of the magnetic transition, is sufficiently precise around M = 0. 
With increasing M, the Coulomb form factor comes into play. In principle, parametrizations of the 
A(1232) Dalitz decay should be checked against the available photoproduction data. The quark 
counting rules constrain the extrapolation to the M ^ region. 

In the present work, we apply the parametrizations of Ref. 



46( 1 ■ which are covariant and 



kinematically complete, i.e., formulated in terms of magnetic, electric, and Coulomb transition 
form factors. The M dependence is based on the extended VMD (eVMD) model and constrained 
by photoproduction and electroproduction data for the form factors, by the transition helicity 
amplitudes of the nucleon resonances, and when available by the p- and w-meson decay branchings. 



The authors of Ref. 



45l | applied the parametrizations of Ernst et al. [351]. The same 



parametrization has also been used for the PLUTO event generator of the HADES group j38l |. 
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However, the work of Ernst et al. fsS^ considered the magnetic form factor with no M dependence, 
which is a crude estimation. Furthermore, the kinematic factors of the A(1232) Dahtz decay in 



Refs. 
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38 



45l ] are incorrect, as discussed in Ref. [90(]. 



Therefore, the interpretation of heavy-ion data still suffers from uncertainties unrelated to the 
complexity of heavy-ion dynamics; i.e., the various parametrizations of the resonance decays used 
as input in the transport models do not agree with each other. 
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FIG. 9: (Color online) Dilepton spectrum in C+C reaction at 2.0^4 GeV as predicted by the vacuum 
NRD+eVMD model and compared with HADES data 'is']. The contribution of the different types of 
sources taken into account in the calculation is explicitly shown. 



2. Collisional broadening 

Let us now turn to the introduction of in-medium effects according to the standard treatments 
and address first Fig. \T0\ in which the HADES data are compared with calculations where the 
possible broadening of the vector meson spectral function in medium is effectively taken into 
account through the introduction of a collisional width Ty^^. We present calculations that use a 
linear parametrization of the type Ty^{pB) = F™^ -|- ps/ po^v^^ipo) to estimate the vector meson 
in-medium width Ty^{pb)- 

In a first approximation, we make no additional assumption concerning the energy dependence 
of the in-medium width: i.e., the same energy dependence is assigned to the collisional width as 



to the vacuum width 36|, |5a]. In particular, the vector meson thresholds are kept the same as the 



vacuum ones, namely, 2m7r and 3m,r for the p and uj meson, respectively. The approximation will be 
investigated below. Figure fTUTa) refers to the assumption F*°*(po) = 200 and T^^{pq) = 60 MeV, 
which reflects the estimates of the CLAS and TAPS experiment for the collisional broadening of the 
p and Lo meson, respectively. Figure fTOT b) refers to the assumption F*°*(po) = 250 , ^]^^ipo) = 125 



MeV. The latter reflects the lower limit estimates emerged from the analysis performed in Ref. [36l |. 
where the values of 300 MeV and 200 — 300 MeV, respectly, for the p and io meson widths at an 
average density of 1.5po were extracted from fits to the DLS data. 
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FIG. 10: (Color online) Dilepton spectrum in C+C collisions at 2.0A GeV for different values of the in- 
medium p and LU widths, (a) r*°t(po) = 200 MeV and r;;;"(po) = 60 MeV. (b) Tf^po) = 250 MeV and 
r^°t(po) = 125 MeV. 

We observe a suppression of the peak with respect to the vacuum case, more pronounced in case 
Fig. 110(b) than in Fig. 110(a). However, in both cases, the experimental data are still overestimated 
around M ~ 0.7 GeV, mainly due to the still significant contribution of the A'^*(1535) resonance. 
The Dalitz decay of this resonance plays a dominant role in the determination of the dilepton 
spectrum in the region around the vector meson peak, due to its strong coupling to the u: meson. 
On the one hand, the HADES data seem to favor a smaller contribution of the A^*(1535) resonance 
in the peak region; on the other hand, however, dilepton production data in p+p collision have 
been well described under the same assumptions for the coupling to the A'^*(1535). This shows 
that the contribution of the A/^*(1535) Dalitz decay, significant in elementary reactions and thus 
in vacuum, is partially reduced in heavy-ion collisions thanks to in-medium effects. We conclude 
that the HADES data suggest a stronger in-medium modification of the uo properties than the one 
taken into account in Fig. [TOl 

Let us now investigate the effect of different choices for the energy dependence of the collisional 
width. In this context, we would like to point out that the mere fact of having, and facing, a certain 
freedom in the choice of an energy dependence of the collisional width shows exemplarily the limits 
that such schematic models carry. Such choices can be based on more or less educated guesses. 
However, if microscopic calculations of in-medium effects are performed, energy dependences are 
fixed from theory, which should be fulfilled as a minimal requirement for a consistent investigation 
of vector meson in-medium properties. Obviously a microscopic calculation of the exact energy 
dependence of the collisional broadening is equivalent to a full model calculation of the in-medium 
spectral function. This well be done later on in this work. 

For the moment, we investigate the consequences of various approximations on a schematic 
level. For this purpose, we extract possible energy dependences of the collisional widths on the 
basis of qualitative considerations and consider the influence on the shape of the final dilepton 
spectrum. Schematically, the collisional broadening that a vector meson acquires is attributed 
to an absorption process of the type V + N^R^'k + N. To simplify, we approximate the 
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corresponding phase space by the phase space for the process M + tti-tv — > 1^1.,^ + rriN and assume 
that the resonance decay proceeds through a p wave. This latter freedom demonstrates again the 
hmits of such schematic procedures. However, since this estimate is only qualitative, let us neglect 
for the moment these refinements. 
One obtains 



-icoU 



{M,p)=T'^\mv,p) 



\M + mN ) 



q{M + mN,'m]\!,mT^ 
q{mv + mN, tun, m-n 



with 



y/[{M + mN)'^ - {rriN + niT^f] [{M + niNY - {rriN - m^) 

q{M + mN,mN, m^j = r 

2 [M + rn-Tv) 



(26) 



(27) 



As one can see, in this approximation, the vector meson threshold is shifted from 2m7r to m-^ for 
the p meson and from SuIt^ to rriTj for the co meson. The choice affects the shape of the lv width 
much more than the shape of the p width. The influence of the choice for the energy dependence 
of the collisional width is illustrated in Fig. [TT]for the case pB = 2pQ and for T^°^{pQ,mp) = 250 
and Tl^^{po,muj) = 125 MeV. In particular, for the to meson, the shift of the threshold leads to a 
large enhancement of the uj width at lower invariant masses. 
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FIG. 11: (a) Imaginary part of the in-medium self-energy — $jSp°'(p, M) = mpT^°^{p, M) of the p meson 
in vacuum (full line) and at p = 2 po for r*°*(po, mp) = 250 MeV (dashed and dashed-dotted lines), (b) 
Imaginary part of the in-medium self-energy — 3Ej|j°*(p, Af ) — mtjr^°*(p, M) of the ui meson in vacuum (full 
line) and at p ~ 2 pa for Tl^^{po,muj) = 125 MeV(dashed and dashed-dotted lines). For both panels, the 
dashed line corresponds to the assumption that the collisional width has the same energy dependence as 
the vacuum width. The dashed-dotted line corresponds to the assumption that the collisional width has the 
energy dependence in Eq. (pS]). 



However, one has to keep in mind that the in-medium p and u widths enter into the expressions 
for the the covariant form factors, see Eqs. (|13p - (jl5p . Their modulus squared determines the width 
r(i? — > N'j*) [Eq. (j24p ]. Thus, only when appreciable differences arise in the covariant form factors 
will the difference in the energy dependence of the in-medium vector meson width be visible in 
the final dilepton spectrum. Now let r*J-^\M) be the p meson in-medium width with an energy 
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dependence analogous to the vacuum width and T*J-^\m) the p meson in-medium width with an 
energy dependence according to Eq. (f26l) . Correspondingly, we set 

(28) 



p[l/2] 
P 



mi 



_p[l/2] 



M2 



M2 



(29) 



We refer now to the uj meson, but the same considerations are valid for the p meson. It can be 
easily realized that \F^\^ « \f}^\^ « 1 when {ml - M'^f > {mj:*^^'^ {i = 1, 2). Thus, the mass 
region where sensible differences between the two cases can be found is typically restricted to the 



=p [21 ^ 

mass region around the vector meson peak. Concerning the p meson, Tp and T'p-''^ are practically 
identical in the region of the vector meson peak, as can be seen from Figllll 



.[1] 



Therefore, we do not expect differences between |-Fp^'p and \F^ 



Concerning the u; meson. 



rj and rj differ substantially in the peak region, although the main differences arise at lower 
masses, i.e., from slightly above m^^ up to slightly above SttItt, because of the different thresholds. 
In addition, one should also consider interference terms of the form Fp F^^ . These terms can, 
however, drive either a constructive or destructive interference, and therefore it is not possible to 
comment on their effect in general within a simple scheme. 

The resulting dilepton spectra obtained for the two choices discussed above are shown in Fig. [T2l 
Here one finds that the contributions from the A resonances, which couple only to the p meson, are 
practically identical in the two cases. Slight differences are visible for the N* resonances around 
the vector meson peak. The differences are more evident in the case of larger values of the widths, 
Fig- [T2r b). However, even in this case, the total spectra differ at most by a factor of 1.3 j96l |. 




0.2 0.4 0.6 0. 

M [GeV] 



0.2 0.4 0.6 0.; 
M [GeV] 



FIG. 12: Dilepton spectrum in C+C collisions at 2. OA GeV for different values of the in-medium p and 
w widths and different choices for the energy dependence of the collision width. The thick lines refer to 
an energy dependence estimated from the V + N^R^n-\-Na,s discussed in the text. The thin lines 
correspond to the same calculations shown in Fig. [TOland are shown for comparison, (a) rp°*(po) = 200 and 
T*-°\pq) = 60 MeV. (b) r*°*(po) = 250 and r;;;''(po) = 125 MeV. 



For a consistent evaluation of the energy dependence resulting from V + N^R^n + N 
processes, one should sum up over all important resonances that couple to the N + V system. 
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each taken with a different weight according to their relative coupling strength, and determine 
for each mode the corresponding angular momentum of the ttN scattering amplitude. Moreover, 
the invariant mass squared of the intermediate resonance would be s = (pjy + p)'^ which leads 
to a dependence on the three-momentum p of the vector meson. It is then clear that such a 
procedure would finally be analogous to the evaluation of the full spectral functions. In fact, the 
V + N^R^TT + N channel is one of the processes consistently included in our calculation of 
the spectral functions, since the A^vr channel is one of the channels entering in the expression of 
the total width of the resonance. 

To conclude, already these first estimates based on the collisional broadening scenario demon- 
strate that the HADES data show clear evidence for a strong in-medium modification of the vector 
meson properties. Figure [T2] demonstrates, on the other hand, that though the two different choices 
under discussion lead to significant deviations of the vector meson widths, particularly concerning 
their threshold behavior, such effects are washed out to a large extent in the final spectra. However, 
the same argument demands a theoretical description that is as precise as possible; i.e., realistic 
spectral functions should be applied rather than pushing schematic models too far. 

3. Dropping mass scenario 

However, before adopting realistic spectral functions, we want to investigate the dropping mass 
scenario d la Brown-Rho, which has been widely used in the literature. Thus, we performed calcu- 
lations for an in-medium scenario that differs from the previous one by the additional assumption 
that the vector meson mass scales with density according to a my = my(l — apB/po) law, with 
a = 0.2. The results are shown in Fig. 1131 where Fig.[T3ja) refers to the choice T^°^{po) = 200 MeV 
and T^fipo) = 60 MeV, and Fig. to the choice r*°*(po) = 250, T^°\po) = 125 MeV. The 
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FIG. 13: (Color online) Dilepton spectrum in C+C collisions at 2.0 AGeV for different values of the in- 
medium p and to widths when an in-medium vector meson mass my = my{l — aps/pa) is introduced, (a) 
r*°*(po) = 200 and r^°'(/0o) = 60 MeV. (b) r*°'(po) = 250 and r^°*(po) = 125 MeV. 

inclusion of a dropping in-medium vector meson mass results in a global shift of the vector me- 
son spectral strength to lower masses. Thus, the corresponding theoretical spectrum is enhanced 
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at lower invariant masses resulting in a sizable overestimation of the experimental data in the 
0.4 < M < 0.7 GeV mass region. At the same time, the experimental data are underestimated 
in the region around and above the vector meson peak because of the lack of spectral strength 
around the (vacuum) vector meson peak. Note that the same underestimation of the vector meson 
peak was observed when adopting the dropping mass scenario to the recent high resolution CERES 
data [91*]. The CERES analysis focused, however, only on the in-medium p meson. 

Concerning the low mass region, m,r < M < 0.4 GeV, the presence of additional strength moves 
the spectrum closer to the experimental data in the mass region M ~ 0.3-0.4 GeV. However, this 
region remains slightly but systematically underestimated. In summary, one can conclude that 
a naive Brown-Rho scaling is too schematic in order to explain the spectrum. This finding is 
consistent with the previous theoretical analysis of the DLS data at lA GeV [3J, |35 ] . 



4- In-medium spectral functions 



Let us now pass to the investigation of in-medium properties based on the in-medium self- 
energies of the vector mesons calculated within NRD+eVMD. First, we present in Fig. [T^a) the 
dilepton spectrum obtained with "first iteration" p and uj spectral functions, i.e., neglecting the 
in-medium modification of the nucleon resonance widths. 
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FIG. 14: (Color online) Dilepton spectrum in C+C collisions at 2.QA GeV resulting from the inclusion of 
p- and w-meson spectral functions calculated within the NRD+eVMD model. The spectral functions affect 
the branching ratios for the Dalitz decays of the baryon resonances, as explained in the text, (a) Inclusion 
of vector meson self-energies determined from vacuum nucleon resonance properties, (b) Inclusion of vector 
meson self-energies calculated in a self-consistent iteration scheme that takes into account the in-medium 
modifications of the nucleon resonance widths induced by the in-medium spectral functions of the vector 
mesons. 



The spectral functions induce a depletion of the theoretical spectrum in the mass region 0.45 < 
M < 0.75 GeV which is not supported by the data. The result can be better understood with the 
help of Fig. [T5t which shows the corresponding p and to contributions which enter into the nucleon 
resonance form factors and determine thus the dilepton production rates. The form factors are 
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determined at saturation density and twice saturation density, pb = Po and pB = '^Po, in both 
cases for a vector meson at rest relative to the nuclear medium (dashed lines). 

The complex structure of the mesonic self-energies is clearly reflected in the form factors, which 
no longer preserve the simple Lorentzian-like shape typical for the vacuum. In particular, for the 
p as well as for the uj we observe a strong minimum around 0.5 ^ M ^ 0.6 GeV between two 
maxima at 0.4 < M < 0.5 GeV and M ~ 0.8 GeV. 

The particular shape of the form factor is determined by the interplay of both the real and 
imaginary part of the self-energy. However, switching off the real part of the self-energy, we 
observe that the depletion of the form factor between M ~ 0.5 and 0.8 GeV is mainly caused by 
large values of the imaginary part of the self-energy in this region. The latter is shown in Fig. [T6l 
This increase is due to the strong coupling to specific resonances, i.e., the A^*(1520) for the p meson 
and the A^*(1535) for the lo meson. The corresponding bump structure is a typical feature for this 
class of models which couple vector mesons to resonance-hole states. 

The inclusion of the in-medium resonance properties caused by the vector meson spectral func- 
tions, i.e., self-consistency, reduces the imaginary part of the self-energy in this region (see Fig. [IB]) . 
In the case of the uj meson, for example, the reduction at M = 0.57 GeV is about a factor of 2.5. 
As a consequence, the form factors, shown in Fig. [TU are enhanced. 

This has an effect on the dilepton spectrum. The spectrum obtained with self-consistent spec- 
tral functions is shown in Fig. [T^bl. The inclusion of the in-medium properties of the nucleon 
resonances moves the theoretical spectrum closer to the experimental data in the mass region 
0.45 < M < 0.75 GeV. This demonstrates the importance of higher order effects, i.e., taking in- 
medium modifications for the nucleon resonances into account when the vector meson properties 
are described by the coupling to nucleon-resonance hole states. 

For the mass region M > 0.4, we conclude that the parameter-free determination of the in- 
medium dilepton spectrum, performed within an approach that attempts to describe simultaneously 
and with the same model parameters 97[ the phenomena of dilepton and vector meson production 
as well as their in-medium modifications, gives a reasonable description of the experimental data. 

However, some data points remain still underestimated. This suggests that the NRD+eVMD 
model predicts a too strong absorption of vector mesons. One possible reason for the present 
underestimation of the experimental data is the use of some poorly constrained eVMD model 
parameters, in particular the RNto couplings. Probably the most relevant case is the A^*(1535) 
resonance, with its strong coupling to the uj meson predicted by the eVMD model though a decay 
of this resonance to Nuj has not been measured yet. Another reason might be that the w-meson 
spectral function in particular is not normalized in the mass region of our interest. The violation 
of normalization ranges from about 30% at p = po to about 45% at p = 2pQ. 9a] In principle, 
this represents no inconsistency, since spectral functions must satisfy the sum rule in the entire 
invariant mass range (up to M = cxd) and not necessarily already in the finite mass interval in 
which we work. This effect should, however, be investigated in future work. 

Regarding the low mass region, ttt-tt < M < 0.4 GeV, the introduction of in-medium spectral 
functions does not provide a solution for the underestimation of the experimental data. On the 
contrary, because of the finite value of the imaginary part of the self-energy at M ~ for high vector 
meson three- momenta p [^^^^(M = 0) 7^ for p 7^ 0], at high momenta we have \ Fy{M = 0)p < 1 
with a consequent reduction of strength. We can therefore conclude that for the explanation of 
the low mass region, one has to take into account additional effects and/or sources. 
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FIG. 15: (Color online) (a) Modulus squared of the p meson contribution to the covariant form factor |Fpp 
at pb — Po (thin lines) and ps = 2po (thick lines), (b) Modulus squared of the u meson contribution to 
the covariant form factor at ps ~ po (thin lines) and = 2po (thick lines). For both panels the 

dashed lines correspond to vector meson self-energies calculated from vacuum nucleon resonance properties. 
The full lines correspond to vector meson self-energies calculated in a self-consistent iteration scheme that 
takes into account the in-medium modifications of the nucleon resonance widths induced by the in-medium 
spectral functions of the vector mesons. Shaded areas indicate the vacuum form factors. 
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FIG. 16: (Color online) Imaginary part of the in-medium self-energy of the (a) p meson and (b) lu meson, 
in vacuum (dashed-double-dotted lines), at = po (thin lines) and ps = 2 po (thick lines). Dashed 
lines correspond to vector meson self-energies calculated from vacuum nucleon resonance properties. Full 
lines correspond to vector meson self-energies calculated in a self-consistent iteration scheme that takes into 
account the in-medium modifications of the nucleon resonance widths induced by the in-medium spectral 
functions of the vector mesons. 



30 



IV. CONCLUSIONS 

We determined the modification of tlie p and uj meson properties in nuciear matter within a 
resonance model and investigated the nonresonant contributions to the vector meson self-energy. 
For both vector mesons, we found a substantial broadening of the width and a significant shift of 
spectral strength down to smaller invariant masses. In particular at small momenta, the coupling 
of the p meson to the iV*(1520)A^~^ state and that of the uj meson to the TV* (1535) A?"""^ state 
lead to pronounced doublc-pcak structures in the spectral functions. In a first approximation, 
the spectral functions were determined from vacuum nucleon resonance properties. Going beyond 
this approximation, the in-medium modification of the nucleon resonance widths induced by the 
modified p and u mesons has been included. This leads to a self-consistent calculation of the vector 
meson spectral functions, which mainly reduces the peaks due to the coupling to A/"*(1520)iV~^ 
and Ar*(1535)iV-i states. 

In a next step, we investigated the influence of different in-medium scenarios for the vector 
mesons on the dilepton production rate in heavy-ion collisions. The dilepton spectrum has been 
calculated exemplarily for the reaction C-l-C at 2.QA GeV for which experimental data have been 
recently released by the HADES Collaboration. 

Already the estimates based on a schematic collisional broadening scenario, i.e., the comparison 
with data, support strong in-medium modification of the vector meson properties. In the dropping 
mass scenario, we found, even when taking additionally into account the collisional broadening 
of the vector meson widths, that the dilepton spectrum overestimates the experimental data at 
invariant masses below the vector meson peak and underestimates them in the region around and 
above the peak. 

Finally, we went beyond the schematic inclusion of in-medium effects and included the vec- 
tor meson properties consistently, i.e., in terms of the in-medium self-energies microscopically 
calculated within our model. We found that self-energies determined from vacuum nucleon res- 
onance properties give a poor description of the experimental data in the invariant mass region 
0.45 < M < 0.75 GeV. On the contrary, the self-consistent iteration scheme provides a reasonable 
description of the data in the same mass region. This demonstrates the importance of consistent 
inclusion of in-medium properties. 

However, for the low mass region (m^r < M < 0.4 GeV) we found that the inclusion of p and 
UJ spectral functions does not improve the theoretical description of the dilepton spectrum and 
experimental data remain slightly underestimated. 

In summary, the investigation represents a first step toward a unified understanding of dilepton 
spectra and vector meson properties in heavy-ion collisions at intermediate energies. The same 
model and the same set of parameters were used to describe the interconnected phenomena of 
dilepton and vector meson production and their in-medium modifications. Forthcoming data, from 
elementary reactions as well as from heavy-ion collisions of heavy systems, will certainly help to 
further reduce still existing model uncertainties. An extension of the present approach to finite 
baryon chemical potential and temperature would further allow one to test the spectral properties 
beyond HADES conditions, e.g., by a comparision with the NA60 data. 
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APPENDIX A: GAUGE INVARIANCE IN EVMD 

The VMD model and its modifications introduce the mixing of a photon with vector mesons 
, (jj, (p, etc. Such a mixing can, in principle, generate finite photon masses and destroy gauge 
invariance. This problem has been solved for the VMD model by Kroll, Lee and Zumomino [93] 
constructing an effective Lagrangian for photons and vector mesons which reproduces the VMD 
predictions. We present first a distinct consistency proof and then show how the method [92] can 
be generalized to the eVMD model. 



1. Final-state interaction (FSI) method 

We start from an effective Lagrangian involving pions interacting with photons. An example 
of such a Lagrangian is the nonlinear a model and, more generally, the chiral perturbation theory 
(ChPT) to a fixed order of the loop expansion. The vector mesons appear as resonances in the 
two-pion scattering channel (p mesons) and the three- meson scattering channel {uj mesons). 



FIG. 17: (Color online) Diagrammatic representation of the FSI of pions (dashed lines) contributing to the 
form factor in the p-meson channel. The photon line is shown as a wavy line. 



Let us consider an absorption of a photon in an isovector channel, as shown in Fig. [T71 Applying 
two-body unitarity and taking into account analyticity (see, e.g., Ref. 93|], Chap. 18), we replace 
the point-like vertex e by ePi{t)/Dj{t) where t = q"^ is the photon momentum squared, Pi{t) is a 
polynomial of the degree /, and Dj{t) is the Jost function defined in terms of the p-wave isovector 
two-pion scattering phase shift 6{t): 



Dj{t) = exp 



5{t')dt' 



TT Jto t'{t-t') 



(Al) 



where to is the two-pion threshold. 

In the no-width approximation, the phase shift accounting for the existence of n resonances is 
given by 



5{t) = Y,^e{t-ml) 



(A2) 



k=l 



32 



where is the mass of the A;th radial excitation of the p meson. Substituting this expression 
into Eq. ()Aip . we obtain 

nt) = m)f{-^. (A3) 
k=i ^ 

The requirement F{t) — > at t ^ oo gives / < n. 

Analytical functions are fixed by their singularities. The representation ()A3p can be rewritten 
in an equivalent additive form 

n 2 

where are some coefficients. The normalization condition F{0) = 1 and the quark counting rules 
impose constraints for c''. 

The effective pion Lagrangian is well defined, since pions are stable particles which exist as 
asymptotic states. In the approach presented above, the problem of gauge invariance does not 
appear, since gauge invariance of the effective Lagrangian ensures a transverse polarization tensor 
of photons and the vanishing photon mass. The vector mesons are resonances accounted for by 
the the final-state interactions. 



2. Effective Lagrangian method 



The vector mesons are unstable particles and do not exist as asymptotic states. Nevertheless, 
the effective Lagrangian method is useful in formulating vector meson effective interactions. Kroll, 
Lee, and Zumomino proposed an effective Lagrangian for the VMD model to illustrate its 
gauge invariance. We extend their arguments for a family of n mesons interacting with photons. 
An effective Lagrangian can be written as 

(A5) 

where is the electromagnetic vector potential, is the kth. p^-meson vector potential, is a 
hadron conserved current, and F^^ = dyA^ — d^j^Ay and G^j^ = dyB^ — dfj_B^. 

Lagrangian (|A5|) is gauge invariant with respect to gauge transformations of the electromagnetic 
vector potential, so the photon interactions with the vector mesons do not violate gauge invariance 
and, in particular, do not generate a photon mass. 

It remains to be shown that the coupling constants and h'' can be chosen such that they 
reproduce the eVMD predictions. For each vector meson, we consider the two diagrams shown on 
Fig. [18l Their sum gives 

n , , 

The spectral functions of the form factors and their asymptotic behavior depend on the type of 
transition. The usual VMD appears for n = 1. It corresponds to asymptotics F{t) ~ 1/t at t ^ oo. 
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If we set =9^, the monopole form factor is reproduced: 



F{t) = (A7) 



For m\ = rup it describes well the pion form factor in the space-like region. The model [9J, |95| of 
the pion form factor, which represents an improvement of the VMD to account for the analyticity 
and two-body unitarity of the pion form factor, and the /j-meson width, works well in both the 
space- and time-like regions. 



FIG. 18: (Color online) Diagrammatic representation of the photon interaction with the electromagnetic 
current using the effective Lagrangian couplings. The first diagram shows the direct photon coupling and 
the second one shows the coupling through the family of n p mesons (double solid line) . The photon line is 
shown as a wavy line. 



In the case of eVMD, we set = c^g^, where are some coefficients. The form factor F{t) 
should decay at infinity, so we obtain 

n 

E = 1- (A8) 

k=l 

Then the usual representation ()A4j) of the eVMD form factors follows. 

The quark counting rules can be satisfied selecting the coefficients . F{t) ~ gives 

n 

J2c'ml = 0. (A9) 

k=l 

Equations (|A8p and (|A9p have a unique solution for n = 2. F(t) ~ requires the existence of 
at least n = 3 vector mesons and so on. 
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